February 8, 2008 23:50 WSPC/INSTRUCTION FILE istanb'lanl 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



O ; PSEUDOSPIN, SPIN, AND COULOMB DIRAC-SYMMETRIES: 

^ ■ DOUBLET STRUCTURE AND SUPERSYMMETRIC PATTERNS 



A. LEVIATAN 

on : Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

ami @phys . huji .ac.il 



- 1—^ 

X 



Received (received date) 
Revised (revised date) 



> 
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Relativistic symmetries of the Dirac Hamiltonian with a mixture of spherically sym- 
' metric Lorentz scalar and vector potentials, are examined from the point of view of 

C ~_ 3 i supersymmetric quantum mechanics. The cases considered include the Coulomb, pseu- 

dospin and spin limits relevant, respectively, to atoms, nuclei and hadrons. 

o 



1. Introduction 

The Dirac equation serves as the basis for the relativistic description of atoms, nu- 
^ | clei and hadrons. In atoms the relevant potentials felt by the electron are Coulombic 

vector potentials. A Dirac Hamiltonian with a Coulomb potential exhibits a fine- 
structure spectrum with characteristic two-fold degeneracy. Relativistic mean fields 
in nuclei generated by meson exchanges 1 , and quark confinement in hadrons 2 nc- 
jJJ ■ cessitate a mixture of Lorentz vector and scalar potentials. Recently symmetries of 

Dirac Hamiltonians with such Lorentz structure have been shown to be relevant 
for explaining the observed degeneracies of certain shell-model orbitals in nuclei 
("pseudospin doublets") 3 , and the absence of quark spin-orbit splitting ("spin 
doublets") 4 , as observed in heavy- light quark mesons. The degenerate doublets as- 
sociated with the relativistic Coulomb, pseudospin, and spin symmetries are shown 
in Table 1. In the current contribution we show 5 that the degeneracy patterns and 
relations between wave functions implied by such relativistic symmetries resemble 



Table 1. Doublet structure in the Coulomb, pseudospin, and spin limits of the Dirac 
equation. The notation (n,£,j) refers, respectively, to the single-fermion radial, 
orbital, and angular momentum quantum numbers of the upper component of the 
Dirac wave function. The Dirac labels are n\ = —(j + 1/2) and K2 = + {j' + 1/2). 
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(n,£,j = e 


+ 1/2) (n - 1,1 + 1, j' = I 


f 1/2) 
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Pseudospin 


(n,i,j = t 


+ 1/2) {n- 1,1 + 2, j' = t 


f 3/2) 


rei + «2 = 1 


Spin 
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2 A. Leviatan 

the patterns found in supersymmetric quantum mechanics (SUSYQM). The feasi- 
bility of such a proposal gains support from the fact that the Dirac Hamiltonian 
with a vector Coulomb potential is known 6 to be supersymmetric. 

2. Dirac Hamiltonian and Supersymmetric Quantum Mechanics 

The essential ingredients of SUSYQM 7 are the supersymmetric Hamiltonian 
H = ( Hl Q and charges Q_ = (°°), Q + = (q^) which S enerate the 

supersymmetry (SUSY) algebra [H,Q±] = {Q±,Q±} = 0, {Q-,Q+} = H. The 
partner Hamiltonians Hi = L^L and Hi = LL^ satisfy an intertwining relation, 

LH X = H 2 L , (1) 

where in one-dimension the transformation operator L — ^ + W(x) is a first- 
order Darboux transformation expressed in terms of a superpotential W{x). The 
intertwining relation ensures that if ^i is an eigenstatc of Hi, then also \&2 = L^!\ 
is an eigenstate of Hi with the same energy, unless L^i vanishes or produces an 
unphysical state {e.g. non-normalizable). Consequently, as shown in Fig. 1(a), the 
SUSY partner Hamiltonians Hi and H2 are isospectral in the sense that their 
spectra consist of pair-wise degenerate levels with a possible non-dcgcncratc single 
state in one sector (when the supersymmetry is exact). The wave functions of 
the degenerate levels are simply related in terms of L and . Such characteristic 
features define a supersymmetric pattern. In what follows we focus the discussion 
on supersymmetric patterns obtained in selected Dirac Hamiltonians. 

The Dirac Hamiltonian, H, for a fermion of mass M moving in external scalar, 
Vs, and vector, Vy, potentials is given by H = a-p + (J(M + Vs) + Vy, where 
a, j3 are the usual Dirac matrices and we have set the units h = c = 1. When 
the potentials are spherically symmetric: Vs — Vs(r), Vy = Vy(r), the operator 
K = — P(cr ■ £ + 1), (with <r the Pauli matrices and £ = —ir x V), commutes 
with H and its non-zero integer eigenvalues k — ±(j + 1/2) are used to label the 
Dirac wave functions m = r~ 1 (G K [Y e x]m,iF K [ Y e > x]m ) ■ Here G K (r) and F K (r) 
are the radial wave functions of the upper and lower components respectively, Ye 
and x are the spherical harmonic and spin function which are coupled to angular 
momentum j with projection m. The labels k = —(j + 1/2) < and £' = £+ 1 hold 
for aligned spin j = 1+1/2 (si/2,P3/2) etc.), while n = (j + 1/2) > and £'=£—! 
hold for unaligned spin j = £ — 1/2 (pi/2,d 3 / 2 , etc.). Denoting the pair of radial 
wave functions by 

•-_(£); « 2 » 

the radial Dirac equations can be cast in Hamiltonian form, 



A(r) = V s + Vy , Z(r) = V s -V v . (3) 
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Fig. 1. Schematic supersymmetric patterns in (a) SUSYQM and in the (b) Coulomb, (c) pseu- 
dospin, (d) spin, limits of the Dirac Hamiltonian. In (a) Hi and H2 have identical spectra with 
an additional level for Hi when SUSY is exact. Spectroscopic notation n£j in (b)-(d) refers to 
quantum numbers of the upper component, and k, N, £ are defined in the text. In (b) the radial 
nodes n are related to n r by n r = n (n r = n + 1) for k < (k > 0), and only the E^ K branch 
is shown 5 . 



In analogy to Eq. (1) we now look for radial Dirac Hamiltonians H Kl and H K2 which 
satisfy an intertwining relation of the form 

LH K1 =H K2 L. (4) 

Following Ref. [8] we consider a matricial Darboux transformation operator 



L = A(r)—+B(r) 



(5) 



where A and B arc 2x2 matrices with r-dependent entries Aij(r), Bij(r). Re- 
lations (4) and (5) should be regarded as a system of equations for the unknown 
operator L and the so- far unspecified potentials in H K (3). The transformation 
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4 A. Leviatan 

operator, when found, connects the two degenerate states 

L$ K1 =C$ K2 (6) 

and imposes relations between their respective components. 

In the usual application of SUSYQM, one starts from a solvable Hamiltonian 
Hi and uses the intertwining relation to obtain a new solvable Hamiltonian H 2 ■ In 
the present case we employ a different strategy, namely, insist that both partner 
Hamiltonians H Kl and H K2 be of the form prescribed in Eq. (3) with the same 
potentials, and look for solutions of Eq. (4) such that the potentials are independent 
of k. We find 5 three physically interesting solutions which require «i + k 2 = 0, 1, — 1 
and lead to the Coulomb, pseudospin, and spin limits respectively. 



3. The Coulomb limit (/ti + k 2 — 0) 

The solution of Eq. (4) with m + k 2 = fix the potentials to be of Coulomb type 

V S = ^ , W = ^, (7) 
r r 

(omitting constant shifts) with arbitrary strengths, as, ay - In terms of the quan- 
ti ties 771 = (as M + a v E)/\, 772 = (a s E + a v M)/\, A = \/M 2 - E 2 , 7 = 
a/k 2 + otg — Qy , the quantization condition reads: 7 + 771 = — n r (n r = 0, 1, 2, . . . ), 
and leads to the eigenvalues 



Etl/M = -*sav±(n r+1 W(n r + ^-a s + ai / _ 
" r,K/ a\ + (n r + 7) 2 v ' 

The K-dependence enters through the factor 7. The spectrum consists of two 
branches denoted by superscript (+) and (— ). The eigenfunctions are 

1/2 




N r(2 7 + i) 



(9) 



2Mn r \7] 2 (K + r]2 

where E = E^ K and Fi = F(-n r ,2-f + l,p), F 2 = F(-n r + l,2j + l,p) are 
confluent hypergeometric functions in the variable p = 2Xr. The states and energies 
in each branch are labeled by (n r ,K). It is also possible to express the results in 
terms of the principal quantum number N defined as N = n r + |/t|, (N = 1,2,...). 
For n r > 1 the eigenvalues in each branch are two-fold degenerate with respect 
to the sign of k, i.e. E^] K = E^} and E^] K = E^} ■ For n r = there is 
only one acceptable state for each k, which has k < for the (+) branch and 
k > for the (— ) branch. Equivalcntly, for a fixed principal quantum number N, 
the allowed values of k are k — ±1, ±2, . . . , ±(N — 1), — N for the (+) branch and 
k = ±1, ±2, . . . , ±(N — 1),+N for the (— ) branch of the spectrum. 

Focusing on the set of states with K\ = — «2 = k, the levels are separated 
according to the value of |«| = j + 1/2. For fixed k, E^ k is an increasing function 
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Dirac- Symmetries: Doublet Structure and Supersymmetric Patterns 5 

of n r and, as shown in Fig. 1(b), for each value of j we have a characteristic 
supersymmetric pattern. There are two towers of energy levels, one for — \k\ (with 
n r = 0, 1, 2, . . .) and one for (with n r = 1,2,.. .). The two towers are identical, 
except that the £^ + i level at the bottom of the — \k\ tower is non-degenerate. 
Similar patterns of pair-wise degenerate levels with ±k appear also in the (— ) 
branch of the spectrum. However, since for fixed k, E nr ] K is a decreasing function of 
n r , the non-degenerate E^ ^ , , level is now at the top of the +\k\ tower, resulting 
in an inverted supersymmetric pattern. The transformation operator is given by 

(d | g + j Ma+ _as_d_,av_ \ 
dr ' r ' Ki Ki dr ' r \ / 1 

as_d__av. d e_ Ma- I ) V W ) 

Ki dr r dr r Ki / 

where e± — K\ + asa±/ n\ and a± = (as ± ay)- Relation (6) is satisfied with 
C = ^-y/n r (7 — 771) and k\ + n 2 = 0. The operator L connects degenerate states 
with (n r > 1,±k), and annihilates the non-degenerate states with n r = 0. The 
condition k± + k 2 = determines the angular momentum couplings in the full 
Dirac wave functions of the doublet 

l(G Kl [Y tX ]W\ _ 1 (G K2 [Y e+lX ]&\ 

* Kl<0 ' m -r{iF K AY i+lX M) ^ >0 ' m -r{iF K2 [Y eX W J ^ 

where m = -(I + 1) < 0, j=£+l/2 and n 2 = +(£+ 1) > 0, f = £+1/2. 

The explicit solvability and observed degeneracies of the relativistic Coulomb 
problem are related to the existence of an additional conserved Hermitian operator 5 

B = -iK l5 (h - /3ikf) + ( a v M + a s H ) , (12) 

which commutes with the full Dirac scalar and vector Coulomb Hamiltonian, H , but 
anticommutes with K. This operator is a generalization of the Johnson-Lippmann 
operator 9 applicable for as = 0. 

4. The pseudospin limit (/ti + /t 2 = 1) 

The solution of Eq. (4) with K\ + k 2 = 1 requires that the sum of scalar and vector 
potentials is a constant 

A(r) = V s (r) + V v (r) = A . (13) 
In this case the transformation operator is given by 







d K2 



Relation (6) is obeyed with Ki + k 2 = 1, C = b(M + A — E), E = E Kl = E K2 and, 
consequently, the radial components satisfy 

dG Kl Ki _ dG K2 k 2 

—\ 1 "Kl — —, 1 «K 2 1 

dr r dr r 

Fki = F K2 . (15) 
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6 A. Leviatan 

The condition K\ + n 2 = 1 determines the form of the full Dirac wave functions of 
the doublet states 

l(G Kl [Y tX ]$\ _ 1 fG K2 [Y l+2X }^\ 

*" 1<0 ' m = r [iF Kl[ Y e+lX M ) *-> '™ ^ r [iF K2 [Y e+lX & ) {W) 

where m = -(£ + 1) < 0, j = £ + 1/2 and k 2 = £ + 2 > 0, / = £ + 3/2. 
From Eqs. (15)-(16) we see that the lower components of the two states in the 
doublet have the same spatial wave function. In particular, they have the same 
orbital angular momentum £ = £+!, and the same number of nodes, n. Eq. (16) 
and a theorem 10 concerning the nodal structure of Dirac bound states ensure that 
the corresponding upper components have quantum numbers (n,l,j = £ + 1/2) 
for ki < and (n — 1,£ + 2,j = £ + 3/2) for k 2 > 0. These are precisely the 
identifying features of pseudospin doublets 11 in nuclei. The latter refer to the 
empirical observation of quasi-degenerate pairs of normal-parity shell-model orbitals 
with such non-rclativistic quantum numbers. The doublet structure is expressed in 
terms of the "pseudo" orbital angular momentum, £ = £+!, and "pseudo" spin, s = 
1/2, which are coupled to j — £±s. Such doublets play a central role in explaining 
features of nuclei 12 , including superdeformation and identical bands. For potentials 
with asymptotic behavior as encountered in nuclei, the Dirac cigenstates for which 
both the upper (G K ) and lower (F K ) components have no nodes, can occur 10 only 
for k < 0, and hence do not have a degenerate partner eigenstate (with k > 0). 
These nodeless Dirac states correspond to the shell-model states with (n — 0, £, j = 
£+ 1/2). For heavy nuclei such states with large j, i.e., 0g 9 / 2 , 0/i n / 2 , 0i 13 / 2 , are 
the "intruder" abnormal-parity states which, indeed, empirically are found not to 
be part of a doublet 12 . Altogether, as shown in Fig. 1(c), the ensemble of Dirac 
states with k\ + n 2 = 1 exhibits a supersymmctric pattern of twin towers with 
pair-wise degenerate pseudospin doublets sharing a common £, and an additional 
non-degenerate nodeless state at the bottom of the Hi < tower. An exception to 
this rule is the tower with k 2 = 1 (pi/ 2 states with £ = 0), which is on its own, 
because states with K\ = do not exist. 

For potentials satisfying the condition of Eq. (13), the Dirac Hamiltonian is 
invariant under an SU(2) algebra, whose generators are 13,14 

MU)- ,17) 

Here s M = are the usual spin operators, s M = UpS^Up and U p = This rela- 
tivistic symmetry has been used 3 to explain the occurrence of pseudospin doublets 
in nuclei. Eqs. (15)-(16) are derived by exploiting the fact that in the symmetry 
limit the Dirac eigenfunctions belong to the spinor representation of the pseudospin 
SU(2) algebra. The relations in Eq. (15) between the radial components of the dou- 
blet wave functions, have been tested in numerous realistic mean field calculations 
in a variety of nuclei, and were found to be obeyed to a good approximation, cs- 
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Dirac- Symmetries: Doublet Structure and Supersymmetric Patterns 7 

pecially for doublets near the Fermi surface 15 < 16 . By the above discussion these 
results confirm the relevance of supersymmetry to single-particle states in nuclei. 

5. The spin limit (ki + n 2 — — 1) 

The solution of Eq. (4) with K\ + n 2 = —1 requires that the difference of the scalar 
and vector potentials is a constant 

X(r)=Vs(r)-V v (r) = Z . (18) 

The transformation operator is given by 



/(2M + S + A) + 

\ fr + *7 

It connects the two doublet states as in Eq. (6) with m + n 2 = — 1, C = —b(E + 
M + E ) and E = E Kl = E K2 . The corresponding radial components then satisfy 



G Kl — G, 



K2 



dF Kl Ki _ dF K2 k 2 

—, r Kl - —, ^k 2 • 

dr r dr r 

The condition m + n 2 = — 1 determines the form of the full Dirac wave functions 
of the two states in the doublet 

* - l ( G >«\ Y 'X$ \ * ( ^[YtxW \ (21] 

*" 1<0 ' m -r[iF Hl [Y» lX ]%) * K > >0 ' m -{iF K2 [Y e _ lX W) (21) 

where m = -{£ + 1), j = £ + 1/2 and k 2 = +£, f = 1-1/2. Using Eq. (20) we 
see that the upper components in Eq. (21) share the same spatial wave function, 
and have quantum numbers (n,£,j = £ + 1/2) for m < and (n,£,j = £ — 1/2) 
for k 2 > 0. As shown in Fig. 1(d), the spectrum consists of towers of states with 
Ki + k 2 = —1, forming pair- wise degenerate spin doublets. In this case, none of the 
towers have a single non-degenerate state and hence, the spectrum corresponds to 
that of a broken SUSY 7 . The tower with m = —1 (si/2 states) is on its own, since 
states with n 2 = do not exist. 

For potentials satisfying condition (18) the Dirac Hamiltonian is invariant under 
another SU(2) algebra, whose generators are obtained from Eq. (17) by interchang- 
ing s u and Sft 13 

Mo°J' (22) 

Eqs. (20)-(21) follow from the fact that the Dirac eigenfunctions in the spin limit 
are spinors with respect to this SU(2) algebra. The spin doublets resulting from 
this relativistic symmetry were argued to be relevant for degeneracies observed in 
heavy- light quark mesons 4 . 
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8 A. Leviatan 

6. Summary 

We have examined three symmetry limits of a Dirac Hamiltonian with spherically- 
symmetric scalar and vector potentials, from a supersymmetric quantum mechanics 
perspective. In the Coulomb limit the potentials are 1 jr but their strengths are oth- 
erwise arbitrary. In the pseudospin or spin limits there are no restrictions on the 
r-dependence of the potentials but there is a constraint on their sum or difference. 
These relativistic symmetries lead to degenerate doublets with quantum numbers 
shown in Table 1, and impose relations between the respective doublet wave func- 
tions. The latter relations are precisely the conditions needed for the fulfillment of 
an intertwining relation, Eq. (4), which is the underlying mechanism of SUSYQM. 
The resulting supersymmetric patterns exhibit sectors of pair-wise degenerate dou- 
blets, with a possible nondegenerate single state in one sector. It is gratifying to note 
that the indicated supersymmetric patterns arc manifested empirically, to a good 
approximation, in physical dynamical systems. While previous studies have focused 
on properties of individual doublets in nuclei and hadrons, it is the grouping of sev- 
eral doublets (and intruder levels in nuclei) into larger multiplets, as discussed in the 
present contribution, which highlights the fingerprints of supersymmetry present in 
these dynamical systems. 
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